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SUMMARY 


The equations of motion for a liquid layer of uniform thickness spread 
over a rigid circular cylinder suhjected to lateral acceleration are derived. 
The liquid is assumed to he inviscid and incompressible and the analysis is 
restricted to the linearized (small deflection) equations of motion. The 
effect of surface tension of the liquid is included. 

A modal solution of the equations of motion is obtained and the stability 
criterion for the liquid layer is derived. In addition, the criterion for 
maximum instability, or droplet formation, is derived. Both of these criteria 
are presented graphically as a plot of axial (along the axis of the cylinder) 
wave length versus Bond number. Correlation between a previous approximate 
theory and the present theory is presented. The available experimental data 
are compared with the theories . 


INTRODUCTION 


The transportation, storage, and utilization of liquids "under near zero 
gravity conditions present many unanswered questions related to the operation 
of spacecraft and space stations . As the effective gravity field is dimin- 
ished, the static and dynamic behavior of a liquid free surface is drastically 
changed since the interface surface tension forces become predominant. 

Most investigators^ to date have considered only the hydrostatic or 
equilibrium configurations of the free surface , while few have endeavored to 
increase the quite limited knowledge of the stability and dynamic character- 
istics of the free surface . The purpose of this report is to investigate in 
some detail a geometrically simple three-dimensional hydrodynamic problem and, 
thereby, to obtain qualitative and quantitative information which may lead to 
a basic understanding of this and more complicated problems . 

The physical problem may be likened to the droplets formed on the under- 
side of a telephone wire dixring a rain storm. The droplets appear to occur 
in a regular distribution along the length of the wire. The same phenomenon 
appears when the coolant tubes of a condenser are located in a horizontal 
position. One might wonder what (if any) change would occur in the droplet 


^Extensive bibliographies are contained in references 1 and 2. 



distribution if the condenser were in a reduced gravity field and^ in fact^ 
if possibly under these conditions no droplets would form. Such a consider- 
ation has significance in connection with condensers carried on space 
vehicles . 

This analytical work is an extension to three space dimensions of that 
reported by Anliker and Beam (ref. l) which presented the two-dimensional 
analysis of the stability of a uniform liquid layer covering a rigid circular 
cylinder. An approximate solution to the three-dimensional problem was 
obtained by Lee (ref. 3 ) and a comparison will be made between his solution 
and the one contained herein. 

In the following analysis^ the liquid is assumed to be incompressible and 
inviscid and the flow is assumed to be irrotational . The density of the 
mediimi external to the liquid free surface is assimied to be negligible. The 
analysis is further restricted to the linearized (small deflection) equations 
of motion^ a rigorous derivation of which is given in reference 1 . The accu- 
racy of the assumption of an inviscid fluid in the investigation of the small 
displacement stability seems to be validated in part by the work of Bellman 
and Pennington (ref. h ) . They have shown that for flat layers^ viscosity 
does not change the stability criterion but merely decreases the rate of 
growth for unstable wavelengths and causes damping for the stable (oscilla- 
tory) wavelengths. The essential discrepancy between the physical problems 
mentioned above and the mathematical problem is that the physical layer will 
not remain precisely at uniform thickness unless the appropriate external 
nonuniform (variation around circumference of layer) pressure distribution is 
applied. The applicability of this theoretical solution to the physical 
problems is then determined by the deviation of the equilibrium configuration 
from a uniformly thick layer. The available experimental data indicate^ 
however^ that the theory is applicable over a wide range of Bond number. 


SYMBOLS 


-^mk 

a 


ao 


B 




g 

Hm 


constant defined by equation (28) 

radius of equilibrium free surface 

radius of rigid cylinder 
pa^g 

Bond number^ — — 

coefficients of Fourier series for antisymmetrical modes 
coefficients of Fourier series for symmetrical modes 
effective gravity (lateral acceleration of rigid cylinder) 
constant defined by equation (23) 
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modified Bessel fimction of the first kind of order m 



Im 

Km 

I>cr 

m,n 

P 

P 

Pe(6) 

Pij 

Q 

R 

r 

S 

T 

t 

X 

X 

P 

Per 

Pmax 

€ 


integers 

modified Bessel function of the second kind of order m 
vavelength corresponding to neutral stability of liquid layer 
integers 

matrix (eq. (^2)) 
press-ore in liquid layer 

equilihri-um pressure distribution applied at free surface 

constant defined by equation (^4) 

matrix (eq* (42)) 

constant defined by equation (43) 

function introduced for separation of variables 

radius of c-urvatiure of free surface 

cylindrical coordinate (fig. l) 

bounding surface 

surface tension coefficient 


time variable 

function introduced for separation of variables 
cylindrical coordinate (fig. l) 
dimensionless wave number 

value of P corresponding to neutral stability of liquid layer 


value of p corresponding to maximum instability of liquid layer 


ratio of cylinder radius to free surface radius , 


ao 

“aT 


© 


displacement of free surface from equilibrium configiaration 
function introduced for separation of variables 
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0 cylindrical coordinate (fig. l) 

Ai axial vave number 

A 2 separation constant 

p liquid density 

frequency parameter (eq. ( 3 I)) 

0 potential function 

^ potential function in terms of space variables 

n potential of body forces 


ANALYSIS 


The mathematical model is formulated by considering the motion of a 
liquid layer of uniform thickness spread over a rigid nonporous circular 
cylinder (fig. l) - The cylindrical coordinates will be denoted by the 
distance along the axis of the cylinder; 0 ^ the angular displacement measured 
from the vertical; and r^ the radial distance from the axis of the cylinder. 
The radius of the rigid cylinder is aQ and the radius of the equilibrium 
free surface of the liquid layer is a. The radial displacement of the free 
surface from its equilibrium position is denoted by The body forces act- 

ing on the fluid (due to lateral acceleration of the cylinder) are normal to 
the axis of the cylinder and parallel to the plane 0 = 0 as indicated in 
figure 1. It is assumed that the liquid wets the rigid cylinder. 

The potential equation for incompressible irrotational flow is^ in 
cylindrical coordinates^ 

V^O(r^x^0^t) = Opp 7 ^ %0 %QC " ^ ( 1 ) 


where 0 is the velocity potential of the fluid,, t the time variable^ and 
a subscript represents partial differentiation with respect to that variable. 
The dynamic boundary condition at the free surface (r = a + is 


P 


- Pe = T 


r=a+^ 



(2) 


where p is the pressure in the liquid layer,, P 0 the external applied 
equilibrium pressure^ T the surface tension (force per unit length) and 

R 2 the principal radii of curvature of the free surface. The kinematic 
boundary condition of a s-urface may be written 



(3) 



■with S = 0 the equation of the hounding surface • At the rigid cylinder 
surface (Sc = r - ao = O) there is no flow through the boundary 



r=ao 


(i<-) 


and at the interface surface (Sj_ = r- a- ^ = 0) there is no flow across the 
sijrface 


DSj_ 

W 



^ % 
d© r^ 



= 0 


which after elimination of nonlinear terms "becomes 


% 



r=a 


(5) 


The pressure^ is given "by the linearized Bernoulli equation (ref. l) 


p = p(^ - ^2 


( 6 ) 


where p denotes the liquid density^ and the potential of the body forces^ 
is given by 


n = gr cos 0 


(7) 


with g representing the effective gravity field coefficient (or the acceler- 
ation of the rigid cylinder) . 

The linearized expressions for the principal radii of curvature Ri and 
Rs are (see ref. l) 

I- = - (8) 


Ra a V ^ ^ 


(9) 


The appropriate su"bstitutions from equations (6), (7) f (8)^ and ( 9 ) into 
equation ( 2 ) lead to 


dt 


- g(a + p cos 9 


- Pf. = T 


^ + 1 Ti _ i - 1 

dx 2 a V a a ^q 2 


( 10 ) 
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which is a combination, of the kinematic and dynamic boundary conditions which 
must be satisfied at the free surface (r = a + ^) or, to the first approxima- 
tion, at the equilibrium free surface (r = a) . From relation (lO) one may 
obtain the pressure distribution which must be applied to the equilibriimi 
free surface to maintain a unifonn liquid layer; namely 

Pe = - ^ - Pga cos e (U) 


Differentiating both sides of equation (lO) with respect to time and 
substituting S^/3t from expression (5) result in 


$ 

PI ^ cos 0 = T| 

ot yr=a 


*^rxx 


+ ^ (5>r + 


-lr=a 


( 12 ) 


The velocity potential ^ must satisfy differential eqmtion (l) and the 
boundary conditions (4) and (12) . With the separation of the time and space 
variables by the substitution 


^>(r,x,0,t) 


±\fat 

e 


3f(r,x,0) 


(13) 


into equations (l), (4), and (12), the stability problem is reduced to the 
eigenvalue problem 


^rr + I ^r + ^ ^00 + ^xx = 0 


(14) 


p (-orif + g¥r cos 0j = T| 
^r=a 




xxr a2 %0r^ 


-'r=a 


(15) 


= 0 


r=ao 


(16) 


where a is the eigenvalue and if the eigenfunction. If 0 is conplex, or 
real and negative, then f increases without limit as time increases and is 
therefore unstable. If 0 is real and positive, then $ is oscillatory 
with time and the amplitude depends on the initial conditions or external 
disturbances; thixs the motion is stable. The stability criterion requires 
that 0 be real and positive. 


Separation of the space variables by the substitution 

= R(r)X(x)@( 0) (17) 

into equation (l4) leads to the following ordinary differential equations 

X" + AaX = 0 
+ Aa© = 0 
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(18) 

(19) 



R" 



( 20 ) 


where Ai and A 2 are separation constants. (The two-dimensional problem con- 
sidered in reference 1 is obtained by setting Ai = 0.) The function 
obtained in this manner will not, however, satisfy the free surface boundary 
condition given by equation ( 15 ) because of the variable coefficient cos d. 
The following procedirre will be used to construct f as a Foirrier series in 
6. Each term of the series will satisfy equation {ik) and the boundary con- 
dition at the rigid cylinder wall, equation ( 16 ) . The coefficients of the 
Fo\rrier series will be chosen so that the interface bovindary condition ( 15 ) 
will be satisfied. In this manner E is chosen of the form 




00 ^ 
m=o 


Im(^i^ ) 


Im (^i^-o^ 
Km (^iS’o) 


Km(Air) 


,iAix 


(cjjj^ cos mS + bm sin m.d) 


( 21 ) 


where Im and Km represent the mth order modified Bessel functions of the 
first and second kind, respectively, and Cm and bm are constants to be 
determined. The separation variable Aa in equation ( 19 ) was replaced by the 
Integer m to Insure physical continuity at 0 = 0, 2it. A prime has been 
used to denote differentiation with respect to the argument of the function. 
The substitution of I from equation (2l) into the free surface boundary 
condition ( 15 ) results in 


Pi - h 


m=o 

00 

+ g COS 0 A i| 

m=o 


Iin(Aia) - Kni(Aia) 

(Aiao) 

00 I 

Im' (Aia) - 


cos(m0) + "bui sin(m0) 


iAiX 


Kjrx (Aiao) 


^321 ( ^ ) 


^ 00 

m=o 


I„'(Aia) - Km' (Aia) 

(Aiao) 


Cm cos(m0) + sin(m0) 
Cm cos(m0) + bm sin(m0) 


.iAiX 




jiAix 


+ 


1 

a2 


00 



Im (Aia) 


m=o 


Im (Aiao) i/T, ^ 

— Kk (Aia) 

Kin v^x^-o) 


Cm cos(m0) 


+ bm sin(m0) 


eiAix 


00 

m=o 


Im (Aia) - 


^m (Aiao) 
Kin *(Ax3-o) 


V (Aia) 


Cm cos(m0) 


+ bm sin(m 0 ) 



( 22 ) 

With the definitions 


_ Im(^i^)Km* (Ai ap) - Im* (^iQ>o)Kin(Aia) 
^ Im* (Aia)Kin* (Aiao) - Im* (Aiao)Kin* (Aia) 


( 23 ) 
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(24) 




Im'(Aia) - Km'(Aia) 

Km (Aiao) 


Cm 


Z 1, 

Dm = •^T" Cm 


-m 


eq.i;iatlon ( 22 ) may "be simplified to 


p -j - a 2 ^ Hme 
m=o 


lAiX 


Cm cos(m 0 ) + "bm sin(m 0 ) 


(25) 


+ g cos e^A.e 
m=o 


iAix 


Cm cos(m 0 ) + bm sin(m 0 ) 


= T ^ > - Ai^e 
m=o 


iAix 


Cm cos(m 0 ) + bm sin(m 0 ) 


00 


Cm cos(m 0 ) + bm siii(m 0 ) 


m=o 


00 




m-o 


Cjn cos(m 0 ) + bm sin( 0 ) 

( 26 ) 


} 


The equations for Cm are obtained by multiplying equation (26) by 
cos(k 0 )d 0 (k = 0 ^ 1 ^ . . .) and integrating from 0 to which produces 


- I Ck + g^Aie^^i^AmkCm 


= T 


- Ai= + ^ (1 - k^) 


m=o 


It ^iAix- 

(27) 




where 


Amk 


pn 

= / cos 0 cos(m 0 )cos(k 0 )d 0 =0 m 7 ^ k ± 1 


( 28 ) 


= ^ m = k ± 1 
4 


Since eqmtion ( 27 ) is to be valid for all values of x, the common factor 
glAix may be omitted 
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r r 00 

o-Hk - ^ ^ + I sAi^Ani,5m = 0 k = 0,1.2, 

m=o 


In a similar manner the equations for hm are obtained by multiplying 
equation (26) by sin(k0)d0 and integrating from 0 to it; thus^ 

r r n “ 

•|- 1 j^- Ax® + ^ (1 - k2)jJ--bi5. + I gAx^Ainkbin =0 k = 1,2, 

m=i 


( 29 ) 


( 30 ) 


which are identical to those for Cm except for the range of summation index- 
The Cm corresponds to symmetrical modes (symmetrical about 0 = O) and the 
bm corresponds to the antisynmetrical modes . 


the following substitutions 


equations (23)^ 

(29), and (30) dimensionless by 

pa^ 

a = — — a 
T 

(31) 

P = Aia 

(32) 

e = f2 
a 

( 33 ) 

T, pa^g 

® “ T 

( 34 ) 


and the final equations for c^ and bj^^ are 


{^5 


Hk(P) - 


P2 + (k2 _ 1) 




A-mk Cm = 0 k = 0^1^2^ 


m=o 


and 


00 

i Hk(P) - p 2 + (k 2 - l)J^bk - I B^Ajnk^ = 0 


k = 1 , 2 , 3 , 


m=x 


( 35 ) 


( 36 ) 


For purposes of clarity equations ( 35 ) and (36) are displayed in the following 
matrix form 
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0 - [p2 + (0 - 

P 

1 )] 

B 

2 


0 


0 


r 

^0 


r 

0 

- B 
2 

P 

- [P2 + (1 

- 1)] 

B 

2 


0 

I 

Cl 


0 

0 


B 

” 2 

5 H2O) 
P 

- [p2 + ik- 

- 1)] 

B 

- 2 


ca 

"1 

0 

0 


0 


B 

2 

5 S3(|3) _ 

[p 2 + ( 9 -l)] • 


C 3 

1 

0 

. 










0 

'w- J 


and 


( 37 ) 


— 






— 


r > 

r 

5 . [p2 + (1 

P 

- 1)] 

B 

2 


0 


0 


bi 


0 

B 

“ 2 

_ H2(P) 
° 3 

- [p2 + {h - 

1)J 

B 

2 


0 


la 

I 

0 

0 


B 

2 

a ^ 3 (P) 

P 

- [p 2 + (9 

- 1)] 

B 

2 

< 

b 3 

■ < 

0 

0 


0 


B 

2 

- H4O) 

- [p2 + (16 - 1)] 


l 4 


0 











1-. 


^ J 


(38) 

There are two sets of linear^ homogeneous^ simultaneous equations for 
the coefficients cm and hm; therefore a nontrivial solution exists only if 
the determinants of the coefficients of equations (37) and (38) are equal to 
zero. The eigenvalue problem for a has been defined 


. [p2 + (0 - 1)] 


554^ - [p2 + (1 _ 1)] 




B 

2 


Ha(P) 


[p2 + (9 _ 1)] 
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( 39 ) 


.4 


for the symmetrical modes, and 


_ Hi(p) 

g 


+ (1 - 1 )] 


0 


0 


B 


2 


[pa + (4 - 1)] 


B 

2 


0 


B 

2 “ 


5 5a(|l . [pa + (9 _ 1 )] 
3 


B 

2 


- _ [-p2 ^ 1)] 

P 


(^ 0 ) 


for the antisymmetrical modes - 


Stability Boxmdary 


The stability boundary will be defined as the bo-undary between regions 
of the p - B plane which correspond to stable and mstable motion of the 
liquid layer. The boimdary will be denoted by Per s-^id the critical wave- 
length by 


2na 

Per 


(la) 


Lq 2- is the minimum length of cylinder for which the liquid layer will be 
unstable. As will be shown in the next section^ Lcr is necessarily (and^ 
in fact^ not usually) the distance between droplets. 

The stability boundary is found in the following manner. Denote the 
frequency determinants (39) sind (4o) by 

iQa - P| =0 (if2) 
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li 



are defined as 


where Q and P are matrices whose elements and 


q. . = + Si 


1 = J 


= 0 


7 ^ j 


Pij = +[ 3^ + (l^ - l)]P 


= 


i = J 
i = J ± 1 


= 0 


^ j ± 1 , i 7^ J 


(^3) 




Since Q and P are syinmetric matrices and Q is positive definite^ then all 
values of a are real (ref. 5 )* The stability criterion that a be real 
and positive may be replaced by the simpler criterion that a be positive . 

The stability boundary is the maximum value of p (minimum value of L) for 
which 0=0. This is accomplished by substituting 0 = 0 into equations (39) 
and (ho) 


p® + (0 - 1) I 

I + (1 

0 I 

0 0 


1) 


0 

B 

2 


f + {h - 1) 


0 

0 

B 


I 3 ^ + (9 - 1 ) 


3 ^ + (1 - 1 ) 


B 

2 

0 


3 ^ + 

B 

2 


0 0 


0 0 

1) I 0 

32 + (9 - 1) I 

i 3 ^ + (16 


1 ) 


= 0 


(^ 5 ) 


= 0 


(46) 
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and evaluating P(b) by assuming values of B and finding the corresponding 
values of p from these new eigenvalue problems. An 8 x 8 (eight rows and 
eight columns) approximation^ was used for the infinite determinants. The 
results of this mmierical exercise^ which was performed on a digital computer^ 
are presented in figure 2. The even mmibered P’s (Pq^ Pa^ • • •) correspond 
to the symetric modes and the odd numbered (pi^ Ps^ . . .) to the antisymmet- 
ric modes. The stability boundary Per is identical to Po since above this 
boundary (P > Po)^ all values of a are positive or^ correspondingly^ P is 
positive definite . 

The symetrical and antisymmetrical mode shapes (0 variation) corre- 
sponding to the lowest ordered p (Pq and Px? respectively) are plotted in 
figure 4. These mode shapes were obtained by substituting the appropriate 
values for P^ B^ and a (a = O) into equations (37) and ( 38 ) to obtain the Cn 

and bn* The 0 variation of the mode shapes was then obtained by the summa- 
r 8 

tions Cn cos n0 and ^ ^ bn sin n0. 

n^o n=i 


An appropriate experimental check on the stability boundary would be a 
systematic increase of cylinder length (keeping the interface radius constant) 
until droplets appeared. Another approach would^ of course^ be to keep the 
length fixed while increasing the interface radius . Experimental data of this 
type apparently have not been obtained. 

Considering thin layers (ao/a ~ l) and assuming a symmetrical mode shape ^ 
Lee (ref. 3 ) obtained an approximate expression for the stability boundary; 
namely^ 

Per ~ (^7) 

Lee’s approximation and the present theory are compared in figure 5* 

B < h the character of the two solutions is quite different. Lee’s solution 
implies that if B < 1^ the layer is stable. The present theory shows Per 
intersecting the p axis at p = 1 which is in agreement with Rayleigh’s 
(ref. 6^ p. 473 ) solution for a liquid jet (no effective gravity^ B = 0^ and 
no rigid cylinder^ ao = O) where the critical wavelength is Ler = 2ata. 


Droplet Formation 


Although the stability boundary^ Per^ separates stable and unstable 
values of p^ it does not necessarily determine the wavelength at which drop- 
lets will form. If the layer is unbounded in the x coordinate (i.e.^ the 

^The convergence of the P’s^ as the size (n X n) of the approximating 
determinant is increased, is illustrated in a plot of p versus l/n for the 
symmetrical modes (fig. 3)* The 8x8 approximation is equivalent to retaining 
a sixteenth-order polynomial for p or an eighth-order polynomial for p^. 

The real roots of the polynomial are the desired P values and are, therefore, 
those values plotted in figure 3 . 
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length of the cylinder is much greater than the radius of the interface 
surface)^ there is a value of to he denoted Pmax^ which a has a 
maximum negative value (i.e.^ a maximum instability).^ Theoretically this is 
the value of P at which droplets will appear.^ The validity of this linear 
analysis in predicting the droplet formation must he ascertained hy a nonlin- 
ear (large deflection) analysis or hy experiment. The experimental work of 
Lee (ref. 3) has shown that linear theory is siafficient to predict droplet 
formation for B > however^ for B < 4 experimental data are^ at present^ 
insufficient for definite conclusions to he drawn. 

The Pmax curve was evaluated numerically for two ratios of rigid 
cylinder radius to interface radius (e = ao/a). The values^ € = 0.0, 0.9^ 
were chosen to he typical of the whole range 0 < e <1. The negative maximum 
value of a was obtained hy taking an 8 X 8 approximation for the frequency 
equations ( 39 ) . The resulting Pmax^ (dd/dp) = 0, curves are shown 

in figure 5* In this same fig-ure Lee's approximate curve is shown for com- 
parison. The analytical expression for Lee's solution is 

Pmax ~ ^ (^8) 

The experimental data presented hy Lee are plotted in figure 5- Fon 
B > theory and experiment agree well,^ however, for B < 4 the experimental 
data are not so complete as to warrant final comparison; therefore, more 
experimental data points are required in this region. 


For the purposes of clarity, a sketch of d (p) for constant B is 
presented in figure 6. The Pej. as well as the Pmax value is indicated. 

The sketch corresponds to ^/B a 3- 

^The argument is that if the response is proportional to a series of 
exponential terms (e^*^ ) and the system is given a general disturbance, then 
the term with the largest exponential factor h will eventually predominate. 

^ee's data was for thin layers. Although he did not specify the thick- 
ness, photographs of the experiments indicate that € = O .9 is a reasonable 
assumption. 

It shoifLd he noted that in his experiments, Lee used a wick at the top of 
the cylinder (0 = O) which provided a boundary condition not taken into 
account in the theory. Conclusive experimental data should he obtained with- 
out obstructions in contact with the cylinder. It should also he remembered 
that the criterion for droplet formation (from the analysis of Lee and the 
present analysis) is based on the assumption of an infinitely long cylinder. 
Thus when the distance between droplets is of the same order of magnitude as 
the length of the cylinder, the end effects may he important. It is interest- 
ing to note that for the two data points corresponding to B ~ 1.2 and 1.6 
(fig. 5 )^ the ratios of test cylinder length to distance between droplets are 
integer values (7 and 8, respectively). In order to meet the physical 
requirements at the ends of the tube, the mmiber of wavelengths must, of 
course, be integer values. 
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CONCLUDING REMAEKS 


The stability criterion for a liquid layer of uniform thickness spread 
over a rigid circular cylinder subjected to lateral acceleration has been 
established according to the linear theory. The criterion for maximum insta- 
bility^ or droplet formation^ of the unbounded (infinitely long) liquid layer 
has been determined. No experimental data are presently available for the 
stability boundary. Experimental data for droplet formation are available for 
Bond numbers 1 < B < 100 . For B > 2 ^ theory and experiment agree well for 
droplet formation. Since the experimental data for B < 2 are quite l im ited 
and somewhat questionable (because of the ratio of the test cylinder length to 
distance between droplets), more experimental data on droplet formation are 
required if the comparison of theory and experiment is to be conclusive. The 
analysis contained in this report is readily adaptable to the study of a 
liquid layer spread over the inside of rigid circular cylinders , such as space 
vehicle propellant tanks exposed to low gravity environments. The stability 
criterion (fig. 2) is directly applicable but the maximum instability crite- 
rion (fig. 5) must be recalciilated. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field., Calif., June 2 , 1964 
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Figure 2.- StaMe region of layer as a fmction of dimensionless cylinder 

length and Bond numher. 
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Figure 3.- Plot of dimensionless wavelength as a function of the size of the 
determinant approximating the infinite determinant in equation (^5) • 
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(a) Symmetrical modes. 



(l) Ant isymmetr leal modes. 

Figure 4.- Symmetrical and ant isymmetr ical mode shapes (6 variation) corre- 
sponding to Po ^or various values of B. All modes are normalized to the 
same constant at their maximum values. 
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”The aeronautical and space activities of the United States shall he 
conducted so as to C07j tribute ... to the expansion of huma^t knowl- 
edge of phefiomena in the atmosphere ajid space. The Adfnmi strati on 
shall provide for the widest practicable and appropriate dissetnination 
of information concerning its activities a?7d the results thereof.” 

— ^National Aeronautics and Space Act of 1958 
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